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Abstract. Basing on main principles of statistical mechanics only, an exact 
virial expansion for path probability distribution of molecular Brownian particle 
in a fluid is derived which connects response of the distribution to perturbations 
of the fluid and statistical correlations of its molecules with Brownian particle. 
The expansion implies that (i) spatial spread of these correlations is finite, 
(ii) this is inconsistent with Gaussian distribution involved by the "molecular 
chaos" hypothesis, and (iii) real path distribution possesses power-law long tails. 
This means that actual Brownian path never can be disjointed into statistically 
independent fragments, even in the Boltzmann-Grad gas, but behaves as if 
Brownian particle's diffusivity undergoes scaleless low-frequency fluctuations. 



PACS numbers: 05.20.Dd, 05.40.Fb, 83.10.Mj 

Keywords : Rigorous results in statistical mechanics, Kinetic theory of gases and 
liquids, Brownian motion 

''God does not play dice'' ( A. Einstein ) 



1. Introduction 

Molecules in gases and liquids, or free electrons and holes in crystals, etc., can be 
treated as "small Brownian particles" (BP) since their thermal motion undoubtedly 
is "not less random" than motion of the R. Brown's pollen suspended in water. 
Therefore, seemingly, the A. Einstein's reasonings [H [2] can be applied to a small 
molecular-size BP too, again producing the diffusion equation for probability density of 
bp's position. The diffusion equation, in its turn, implies that probability distribution, 
Vo{t, AR) , of displacement of BP, AR , during time interval (0, t) , at long enough t 
tends to the Gaussian distribution: 

Vo{t,AR) VG(t,AR) = (ATTDt)-^/^ cxp {-ATi^ /4Dt) (1) 

Formally, Einstein in [T] assumed that AR consists of many increments which are 
statistically independent in the sense of the probability theory. Thus, from statistical 
point of view, his result is equivalent to the " Zaw of large numbers" discovered by 
J.Bernoulli almost two centuries earlier [3]. 
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The idea of statistical independency, in the form of the "Stosszahlansatz" , or 
"molecular chaos hypothesis" , was also taken by L. Boltzmann as a principle of his 
molecular-kinetic theory of gases [1]. Later in [5J N. Bogolyubov imparted it, after non- 
principal modification, to a theory based on the Bogolyubov-Born-Green-Kirkwood- 
Yvon (BBGKY) hierarchy of equations [3 [H [71 [H [9] . It was violence against own 
brainchild, because the BBGKY equations are quite all-sufficient ones and do not 
need in any add-on (except natural initial conditions, of course). Unfortunately, the 
violence lasts up to now: various (higher-order, non-local, etc.) generalizations of the 
Boltzmann kinetic equation to relatively dense gases [SI El [HI IS] all the same rest upon 
one or another variant of the "Stosszahlansatz" . Although an attempt to break this 
tradition was made in [10] (see also [TT]). 

At same time of the forties, N. Krylov in his book [12j (first published in 1950 
in Russian) argued that probability-theoretic concepts under use in modern kinetics 
generally are incompatible with principles of statistical mechanics. Especially he 
emphasized fallacy of common opinions that ^^probabilities do exist regardless of a 
theoretical construct and full-scale experiments" and that obviously independent 
phenomena should possess independent probability distributions" (italics mark brief 
citations from [12]). According to Krylov, "physical independency" of events in reality 
does not mean their "statistical independency" in theory. 

At present, survivability of the prejudices disclosed by Krylov is the only excuse 
of the "Stosszahlansatz" or similar conjectures. They gave rise to conviction that in 
the limit of "Boltzmann-Grad gas" ("infinitely dilute gas") infinite BBGKY hierarchy 
reduces to the single Boltzmann equation, i.e. the latter presents exact gas kinetics. 
If such was the case, then random walk of a molecular BP (e.g. test or marked gas 
atom) would be made of many statistically independent events and pieces. Then the 
"law of large numbers^'' is in effect, and probability distribution of the Brownian path 
AR in thermodynamically equilibrium gas should have Gaussian asymptotic ([1]). 

But this is not true! In fact, as we will show, Vb(i, AR) has essentially non- 
Gaussian asymptotic possessing power-law tails at AR^/4Z?t > 1 instead of the 
exponential ones, even in equilibrium gas under the Boltzmann-Grad limit (BGL). 

Hence, N. Krylov was right, and statistical independency of colliding molecules or 
statistical independency of different pieces of Brownian trajectory, etc., like statistical 
independency of mathematical dice tosses, exists in imagination only but not in 
physical reality. If reinterpreting the well-known words said by Einstein, one can 
say that he was playing dice in [1] but God does not play dice. 

A substantiation of these statements below (see also [T31 dH [H]) will be done 
"at very thermodynamical level" basing on only determinism and reversibility of 
Hamiltonian microscopic dynamics and besides general notions about many-particle 
distribution functions and correlation functions [3 [7l [8] and main fluid properties. 

In spite of so abstract approach, remarkably, our conclusions will be in full 
qualitative agreement with result obtained in [16j by means of crucial approximation 
of the BBGKY hierarchy and then its direct solving under BGL: 

Vo{t, AR) ^ {AwDty 3/2 r(7/2) (1 + AR^ ADt)- (2) 

(in [in] designation Wi was used in place of Vb). The diffusivity D here, as well 
as in ([1]), is defined by J AR^ Vb(t, AR) dAR — > 6Dt , while the arrow everywhere 
denotes asymptotic at t much greater than BP's mean free-flight time. 

Our plan is as follows. In the beginning, a kind of virial expansion for Vo{t, AR) 
is derived. It connects, from one hand, first- and higher-order derivatives of Vo{t, AR) 
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with respect to gas density and, from the other hand, pair and many-particle joint 
correlation functions (CF) of BP and gas. The CF in their turn are definitely related 
to usual distribution functions (DF). Since the latter by their sense are non- negative, 
the virial expansion results in a series of differential inequalities to be satisfied by 
Vo{t, AR) . The first of them, eventually, restricts a steepness of Vb(i, AR)'s tails and 
clearly forbids their exponentially fast decrease. At the end of the paper, its relation 
to the problem of 1/f noise is commented. 

2. Virial expansion of Brownian path probability distribution 

2.1. Full-scale experiments and fluctuation-dissipation relations 

We will consider a system consisting of ^ 1 atoms in volume f2 plus one more 
corpuscular "Brownian particle" (BP), under the thermodynamical limit: N oo , 
J7 ^ (X) , N/n = vq — const. Initially, at time t — position of BP R(t) is 
definitely known. We are interested in the already mentioned distribution Vb(t, AR) 
of consequent BP's path AR ~ R(i) ~ 1^(0) ■ Especially (see Introduction), in 
thermodynamically equilibrium Boltzmann-Grad gas where a fortiori any "collective" 
or "hydrodynamic" contributions to chaotic motion of particles disappear. At the same 
time, as far as possible, we want to take in mind also finite-density and non-equilibrium 
gas and even liquid. 

Instead of "art of decomposition" of BP's path into some constituent parts and 
"art of conjecturing" about their "probabilities" , we want to follow N. Krylov (see 
Introduction) and consider a set of mental "full-scale experiments" to see how BP's 
path as a whole is influenced by artificial controllable perturbations of both gas and BP. 
Thus we involve into consideration statistical correlations between BP's path and gas 
atoms. At that, the only "probability" to be specified (at our own choosing, without 
any conjectures) is initial probabilistic measure in the space of states of the system. 
To find other "probabilities" , one should analyze a flow of this measure according to 
the Liouville equation or equivalent BBGKY equations. But a great deal can be found 
from time reversibility of this flow only which is expressed e.g. by the "generalized 
fluctuation-dissipation relations" (FDR) pJl[T8l[T9] . 

Let, firstly, q = {R, ri, ... , rjv} and p = {P, pi, ... , pjv} are (canonical) 
coordinates and momentums of our system, and -ff (q, p) its Hamiltonian in 
absence of its perturbations, so that iJ(q, — p) = iJ(q, p). Secondly, the initial 
probabilistic measure, /3in(q, p) , is the equilibrium canonical one corresponding to 
this Hamiltonian, that is Pin(q, p) = Peq{^.^v) ^ ^xp [— i?(q, p)/r] (omitting a 
normalization factor). Thirdly, suppose possibility that BP is perturbed by a constant 
external force f which is sharply switched on at t = , thus at t > changing the 
Hamiltonian to H{t, q, p) = (q, p) — f • R . 

Then we can either directly write according to works [17] or |19j or derive by their 
methods, like in [T5], the very particular but useful enough FDR as follows: 

(A(q(t),p(0)i?(q(0),p(0))e-^W/^)o = 

= (B(q(i),-p(t))A(q(0),-p(0)))o 

Here £{t) = f • [R(<) — R-(O)] is work made by the external force during time 
interval {0,t) , A(q, p) and B{q,p) are "arbitrary functions", the angle brackets 
(...)o designate averaging over statistical ensemble of phase trajectories of the system 
corresponding to the equilibrium Gibbs ensemble of their initial conditions, and T 
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is initial temperature of the system (or, to be precise, of the ensemble). At any 
concrete trajectory, of course, q(i) and p{t) represent solutions to Hamilton equations 
corresponding to the Hamiltonian 7J(i,q, p) = _ff(q, p) — f • R and hence are some 
functions of the force f . The equality ^ holds also if BP has internal degrees of 
freedom. In absence of external force, when £{t) = , it is so obvious that even does 
not need in a proof. 

Further it is sufficient to confine ourselves by such particular choice as 

^(q,p) = <5(R-R') , 

B(q,p) = r!5(R-Ro)5(P-Po)exp[-Ef=if^(r„P,)/r] (4) 
= <5(R-Ro)<5(P-Po)nf=i[l + ^(r„P,)] , 
where 'ip{r,p) = exp [— [/(r, p)/r] — 1 . 

2.2. Generating functional of distribution functions 

Under choice ^ in the thermodynamical limit right-hand side of ^ takes form 

( i?(q(t), -p{t)) A(q(0), -p(0)) >o = T{t, Ro, -Po, |R'} (5) 

where 

.F{t,Ro,Po, iZ-IR'} = K)(i,Ro,Po|R') + 

oo ^„ „„ n 

+ F„(t,Ro,ri...r„,Po,pi...p„|R')n^(^J'-Pj) (6) 

n=l "''■XP j=l 

Here symbol /J'^p denotes integration over n coordinates ri...r„ and momentums 
Pi ...pn ; function Vb(t, Ro, Po|R') is conditional probability density of finding BP 
at t > at point Rq with momentum Pq under condition that BP had started 
at t = from point R' , and Rg, ri ... r„, Pq, Pi ... p„|R') is joint conditional 

probability density of this event and simultaneously finding some atoms at points Vj 
with momentums Pj under the same condition. 

In respect to atoms, all the Fn are complete analogues of standardly defined non- 
normalized many-particle DF of infinite gas [5]. Instead of normalization, "asymptotic 
uncoupling" of inter-particle correlations takes place: 

F„(...rfc...pfe...) ~> F„_i(...rfc_i,rfe+i... pfe_i,pA;+i...)G,„(pfc) , 
Fi(i,Ro,Po,ri,pi|R') ^ K)(<,Ro,Po|R')G„,(pi) , 

when Ffe — > oo and ri ^ oo , respectively, with Gm (p) being equilibrium Maxwellian 
distribution of atomic momentum, Gm{p) = {2T:Tm)~^^'^ exp {—p^ /2Tm) , and m 
atomic mass. But, because of initial localization of BP, in respect to BP's variables 
all the DF are normalized in literal sense. In particular, 

JVo{t,Ro,PomdPo = yo(i,Ro-R') , 
JVoit,Ro-K')dIio = 1 

In respect to ip , expression J^{i,Ro,Po, V' |R'} represents generating functional for 
these DF quite similar to the functional originally introduced by Bogolyubov [S]. 

By definition of the average (...)o, all the DF represent BP in initially 
thermodynamically equilibrium gas (or, to be precise, in equilibrium Gibbs ensemble 
of identical systems "gas plus BP" ) . Correspondingly, initial conditions to them are 

yo(0,Ro,Po|R') = <5(Ro-R')Gm(Po) , (7) 
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F„(0,Ro,ri ...r„,Po,pi ...p„|R') = 



(8) 



= Fr'(ri...r„|Ro)^(Ro-R')GM(Po)n,"=iG™(p,) , 
^{0, Ro,Po, V|R'} = <5(Ro-R')Gm(Po)^'^'?H0|Ro} , 



where 



0(r) = /V'(r,-p)G„i(p)rfp = /V'(r,p)G„i(p)dp , 



(9) 




°o ,,n rn 



^(^^){0|Ro} ^ 1 + ^ ^ / F(^^)(ri...r„|Ro) [] '^(r.) 



functions (ri... r„|Ro) are conditional equilibrium DF of gas under fixed position 
of BP, T^'^'^^ {4>\Ilo} is their generating functional, and M means BP's mass. 

Of course, all the DF possess translational invariance. Notice also that, firstly, 
ratios Fn/Vo represent conditional DF of gas under condition that both initial position 
of BP and its current state are known. Secondly, al least under condition 



the infinite series are converging and hence the functionals are well defined. 
2.3. BP-gas correlations 

In absence of the external force ( f = ) just listed DF describe Brownian motion in gas 
which all the time remains in exact thermodynamical equilibrium. Nevertheless, at 
t > not only Fn but conditional DF Fn/Vo too are constantly changing along 
with Vq . In principle, their joint evolution is unambiguously prescribed by the 
BBGKY equations together with initial conditions Q and (0) [M HH [TH [Tg] . At 
that, differences 



i^„(i,Ro,ri...r„,Po,Pi...p„|R') - R-o, Po|R') -P^n (ri- r„|Ro) Jlj Gm(p,) 



reflect additional specific statistical correlations between BP and gas which arise just 
due to evolution of Vo(t, Rq, Po|R') ■ Roughly, an origin of their specificity is that they 
are correlations of a current gas state with previous BP's path accumulated during 
all the time interval (0, t) . By this reason we can characterize them as "historical 
correlations" . 

By tradition, any additions to equilibrium (or quasi-equilibrium) DF are 
termed "correlation functions" (CF) [Zl [3 [!]■ We will apply this term also 
to functions what describe the "historical correlations". Let us designate CF 
as I4,(i, Ro, Ti ... r„, Pq, Pi ... Pn|R') and introduce them through their generating 
functional: 



/ |0(r)| dr < oo 



(10) 



.F{t,Ro,Po, ^|R'} = .F(^«){</)|Ro} V{t,Ro,Po, ^|R'} , 
where = 4>{r) is expressed through ip = ip{r, p) by formula and 

V{t,Ro,Po, ^|R'} = yo(t,Ro,Po|R') + 



(11) 



(12) 
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where function Vi represents pair "historical" correlation between BP and atoms. 

In view of the asymptotic decoupling of inter-particle correlations, it is clear that 
all the CF disappear in initial equilibrium state, V^i(0 ,...) = at n > , as well as 
far from BP, i.e. , ... r^... ) ^ at ^ oo . 

2.4- Main relation between correlation functions and response of Brownian path to 
gas perturbations 

Now let us introduce 5'(q,p)= (5(R — Rq) Hj [1 + ^'(I'j: Pj) ] and, under the choice 
(|H), consider left side of ([3]) rewriting it as 



with $ being arbitrary functional of the system's phase trajectory. 

Evidently, (...) in fact represents averaging over new statistical ensemble of initial 
conditions, with new probabilistic measure 



Pi„(q,p) cx pe9(q,p)*(q,p) oc S(R-Ro) exp[-H{q,p)/T -Y.jU{qj,Pj)/T] 



(again normalizing coefficients are omitted). Formally this measure can be viewed 
as also a canonical equilibrium one but in presence of generalized (momenta- 
dependent) external potential t/(q, p) = —T In [1 -I- p) ] . In fact, of course, this 
is thermodynamically non-equilibrium measure since system's Hamiltonian does not 
include such potential. Hence, second right-hand average in (|14p is nothing but 



where V{t,'R.'\ijj,Ilo,Po} is conditional probability density of finding BP at point R' 
under conditions that initially, at t = , it was located at point Rq with momentum 
Po while gas was in such non-equilibrium spatially-nonuniform state what would be 
equilibrium under external potential U (q, p) . Noticing, besides, that 



Equivalently, instead of the artificial external potential U (q, p) , the non- 
equilibrium ensemble which has arisen can be characterized by corresponding 
conditional mean densities of atoms in the /i-space and coordinate space. 




(<5(R(i)-R')^(P(0)-Po)) = V{t,Ryj,Ro,Po}GM{Po) , 



(15) 




(16) 



^{r,p|V',Ro} = (^(5(rj -r)^(pj 



i 



P)> 



i^{r\ (/>,Ro} Gm(p) 



1 + '0(r,p) 
1 + 0(r) 



•{r|<^,Ro} ^ iY^Sir.^r)) 



[1 + 0(r)] 



JlnJ^("'?){0|Ro} 
(50(r) 



(17) 
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At that, /^{r,p| ■(/'jRol/^o has the sense of initial one-particle DF of atoms. 

At last, combining formulas pT|) . and (fTB|) . we come to formally 

exact relation 

T/{t,R'|V^,Ro,Po} Gm(Po) e-f-[^'-^°l/^ = (18) 
= Voit, Ro, -Po|R') + ^0' -Po' Pi -P"!!^') n VXr,, -Pj) 

n=l "''■XP j=l 

It connects, from one (left) hand, probability distribution of BP's path in initially non- 
equilibrium nonuniform gas and, from the other (right) hand, probability distribution 
of bp's path, along with generating functional of correlations between this previously 
accumulated path and current BP's environment, in initially equilibrium uniform 
gas. In case f = , therefore, right-hand side of represents wholly equilibrium 
Brownian motion. 

In other words, relation pS]) connects two sorts of "full-scale experiments": one 
on susceptibility of Brownian motion to perturbations of medium where it takes place, 
and another on its correlations with thermal fluctuations in the medium. In such sense, 
(|18p is typical generalized FDR or "generalized Onsager relation" . 

At V(r,p) = 0, clearly, ^TE\i turns into 

yo(i,R'|Ro,Po) Gm(Po) e^f-[^'-^«l/T ^ i/o(t, R^, _Po|R') , (19) 

where Vb(t, R'|Ro, Po} is density of probability to find BP (in the same gas as on 
the left) at point R' under condition that it started from point Rq with momentum 
Pq . Integration over this momentum yields the classical FDR [18l [191 HH [22] 

yo(t,AR) e-f-^^/^ = Voit, -AH) , (20) 

where AR = R ' — Rq and 

l/o(i,Ro-R') = Jvoit,Ro,'Po\^')dPo 

Of course, all the DF and CF are dependent on the mean gas density Uq (and on 
the force f if any) but for brevity in ([TH]) , as well as before it and almost everywhere 
below, corresponding arguments are omitted. 



2.5. Quasi-uniform gas perturbations, spatial correlations and virial expansion 

Let, firstly, the gas perturbation does not change velocity distribution of atoms, that 
is represents pure density perturbation, ^p{r,p) = (^(r) . Secondly, 0(r) = = const 
inside some sphere |r — Ro| < ^ and vanishes outside it in some suitable way (it 
should be underlined that nothing impedes choosing perturbations to be correlated 
with the point Rq ). Since, according to (fTO)) . 0(r) is absolutely integrable, radius 
must be finite. But it can be as large as we want. At that, factual perturbation of gas 
equilibrium initially is located at |r — Ro| > £, . 

For example, we may take ^ = k v^to , where Vs is speed of sound in our gas (or 
liquid) , to is maximal duration of our "full-scale experiments" , and k > 2 . Then, 
if the external force is not too strong (so that velocity of BP's drift induced by the 
external force is small as compared with Vg), we can be sure that at t < to the 
Brownian motion mentioned in left part of psp practically takes place in equilibrium 
uniform gas with a constant mean density v = const what corresponds to </> = const . 
Indeed, under mentioned conditions |R(t) — R(0)| < Vgt while radial approach of 
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inner front of the gas perturbation (let even Mach front) to the point R(0) — Rq is 
not greater than Vgt . Hence, undoubtedly BP does not feel the front and moves as it 
was in spatially uniform (and thus equilibrium) media at least till t < Iq . This can 
be named "quasi-uniform perturbation" of gas. 

Now, at t < to , in fact both parts of p8)) describe Brownian motions in (initially) 
uniform equilibrium gases but with different values of density. On the right it is vo 
while on the left it equals to i^{r|0, Rq} taken far from BP, at |r — Ro| ^ rs , with 
rs being radius of pair interaction between BP and atoms. We will denote this value 
simply as ly . According to (fT7|) , it is definite function of cj) and the seed density i^o : 

v^v{vo,<f) = i.o(l + 0){l + ;^o0/[Fi''\r)-l]+ (21) 

+ ^ ^ Vi^'Hri,r2,0) - ^^i^«)(ri) - F(^«^(r2) - Fi^'^ri - r^) + 2] + ... } 

Here ^^"'^(r) and Fg^'^^ri, ra, rg) are standard pair and triple DF of equilibrium gas 
with density Vq (i.e. functions what follow from Fi^'^-'(ri...r„|Ro) at Rq^oo). 

Correspondingly, the functional R'|(/), Rq, Pq} simplifies to mere function, 
and, integrating left side of (fT8|) over Pq , for t < to we can write 

/F{t,R'|0,Ro,Po}GM(Po)dPo = K)(i,R'-Ro; ^^) , 

where Vq has exactly the same sense as Vo on the right-hand side of (after its 
integration over Pq ) , and we introduced the density argument v , so that 

Voit, Ro - R' ; = J^o ) = Vo{t, Ro - R') 

On right-hand side of (|18p . under the formulated conditions, a fortiori none 
correlations between BP and gas might propagate out of the sphere |r — Rq] < <^ . 
Therefore in all the integrals (j){r) can be replaced by the constant. Consequently, 
relation ([TH]), after its integration over Po , transforms into 

^o(t,AR; Ki^o,V')) e-'-^^/^ = Fo(t,-AR; i^o) + V ^K(t,-AR; i^o) (22) 

n — 1 

functions Ki at n > defined by 

Vn{t,Ilo — R' ; i^o) ^ '^o / Ki(i, R07 Ti... r„, Pq, pi... p„ I R') dPo (23) 

Jrxp J 

Combining ([22|) with FDR (I20p (which, of course, is valid at any density if equal on 
both sides), we obtain 

^ In 

Vo{t,AR; Hiyo,ij)) - K)(i,AR; i^o) + V ^K(t,AR; i^o) , (24) 

n — l 

where now AR = Rq — R' . 

Relation ((24|) thoroughly may be qualified as "virial expansion" of the BP's path 
probability distribution, Vo(t, AR; v) , similarly to well-known virial expansions of 
thermodynamical quantities [24] or kinetic coefficients [9j. But there is significant 
mathematical difference between them: the two latter expand over absolute value of 
density v while the former in fact over its relative value v jvo ■ 
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2.6. First-order relations 

First-order, in respect to (p , terms of produce 
dVo{t,AK; vo 



Fi(t,AR; i^o) 



rxp 



Fi(t,Ro,Po,ri,Pi|R')dPo (25) 



Here AR = Ro - R' and, according to ^7} and/or ((2T 



dv{vQ,i>) 



dtp 



= V(i + VQ 



(26) 



Here, of course, F!f''^\T) is some function of . The last equality in ([^5]) . with P 
standing for the pressure and the bracket representing compressibility, is well known 
from statistical thermodynamics [24] , 

Notice that the function Vi{t, AR; vq) can be interpreted as total pair correlation 
of BP with whole gas but counted per (elementary volume what falls at) one gas 
atom. Similar interpretation is applicable to higher-order integral correlations (|23p . 

Returning to our basic relation (I18|) in the full phase space, in the first order with 
respect to ip we have 

■<5y{t,R'|^,Ro,Po} 



JV(r,p) 



(27) 



= zyoyi(i,Ro,r,-Po,-p|R') 



With the help of ([T7]), variational derivative with respect to ip can be replaced by 
that with respect to density of atoms in /i-space: 



5V 



(5V(r,p) 



5V 



^=0 <5/i(r,p) 
G™(p')[^^i^^^(r,r'|Ro)-F('^^'(r|Ro)F(^^'(r'|Ro)] 



(28) 



(5/x(r',p') 



with the right-hand variational derivatives taken at /i(r, p) = j/oGm(p)^i''^'''' (r|Ro) , 
i.e. at equilibrium gas with the seed density uo . This generalizes (I26p to arbitrary 
non-uniform perturbations. 

The two latter formulas establish quite rigorous connection of, from one hand, 
pair correlation between previous path of BP and current state of the medium where it 
is walking, and, from the other hand, linear susceptibility of the BP's path probability 
distribution to weak perturbations of the medium. 



3. Ranges of spatial statistical correlations between medium and 
Brownian particle and restrictions on its path probability distribution 

Our previous results clearly show that, firstly, historical correlations between BP and 
medium (gas or liquid) certainly exist, i.e. are not zeros. Secondly, their total values 
(integrated over all momenta and relative distances between atoms and BP as in ((23)) ) 
quite definitely reflect sensitivity of the probabilistic law of Brownian motion to change 
of state of the medium, first of all, to density of its atoms. 

Indeed, let the medium represents three-dimensional weakly non- ideal ("dilute") 
gas in equilibrium (f = 0). Then, undoubtedly, an increase of gas density must 
lead to constriction of the distribution Vb(i, AR;;^). Thus its density derivative 
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9Vo(<, AR; vo)/dvQ is positive at sufRciently small |AR| but negative at relatively 
large |AR| . Since, undoubtedly again, auto-correlation of BP's velocity is integrable 
and therefore BP's chaotic walk is organized as "diffusion" considered already in [Tl[5], 
i.e. AR2(i) cx Dt , / AR^ Vb(t, AR; i^o) ^AR = &Dt , a bound what separates 
"small" and "large" values of |AR| is of order of . Hence, we can state that 

9Vo(t, AR; i/o)/di'o < when z = AR^/4Z?t is significantly greater then unit. 

This implies, according to (pS)) . that total pair BP-gas correlation, Vi {t, AR ; i/q) , 
also is negative at large z, and thus Vi(i, Rq, Po, fi, Pi|R-') somewhere is negative. 
But it can not be "too negative" , because its contribution to right side of (fT3|) 
should not make its left side (which represents probability distribution) negative. 
This means, obviously, that the integral Vi(t, AR; i/q) and thus the derivative 
Vq dVo{t, AH; i'o)/di'Q is bounded below by some negative whose absolute value is 
proportional to Vo(t, AR;i/o) • Such restriction of the derivative, in its turn, means 
that tails of Vo{t, AH; vq) at z ^ 1 can not decrease in too fast way. Anyhow, 
possibility of exponential decrease of Vo{t, AR; i^q) at z ^ 1 and hence the Gaussian 
asymptotic ([T]) become under question. Let us consider this suspicion carefully. 

3.1. Restriction on the BP's path distribution tails. Estimate 1 

At any t , R ' , Rq , Pq and p , let 

/i(t,Ro - R',Po,p) = min t/i(i,Ro,r,Po,p|R') 

r 

Assume that at given i, R', Rq , Pq and p integral of Vi(t, Rq, r, Pq, p|R') over 
r is negative (i.e. Vi(t, AR; uq) < 0). This means that /i(i, Rq — R',Po,p) also is 
negative, and therefore we can introduce effective volume occupied by negative pair 
correlation, ri„eg(i, AR, Pq, p) , by means of 

_ / T/i (^,Ro,Po,r,p|R O dv 
/i(i,AR,Po,p) 

Clearly, this is lower boundary of volumes what could be reasonably attributed to the 
negative correlation. If, in opposite, the integral over r is positive, in the same sense 
it is reasonable to assign ^neg{t, AR, Pq, p) = . Then in any case we can write 

/ yi(t,Ro,Po,r,p|R')dr > a„eg(i, AR, Pq, p) /i(i, AR, Pq, p) (30) 
Next, pay our attention to identity IT^l) . Since the pair DP Fi represents a 
probability density, it must be non-negative, Fi(t, Rq, r, Pq, p|R') > 0. Hence, 

/i(t,AR,Po,p) > -yo(t,Ro,Po|R') G™(p) maxi^i^""^ (31) 

with max i^;^''^'''' denoting maximum of all values of F^'-'^'^' (r|Ro) at those points r 
where (i, Rq, r, Pq, p|R') takes its minimum value ft,(i, Rq — R',Po,p). This 
inequality will be even more so valid if replace max Fi""^^ by absolute maximum 
of Fi : max f['"^'^ maxr f['"^'> (r|Ro) . 

It is useful to notice that Fi(t,Ro,r = Ro,Po,p|R') = 0, i^^^'"'^ (Ro|Ro) = 
and therefore Vi(t, Ro, r — Ro, Pq, p|R') = , if BP and an atom can not be located 
at same point. Consequently, always /i(t, AR, Pq, p) < 0. 

Combining inequalities ([50)1 and PT|) and relation (pS)) . we come to inequality 

vq > ~vq maxF{ ' x (32) 



a„e,(i,AR,Po,p) ^ ■> \u It.\ L (29) 



j j a,eg(<,AR,Po,p)Fo(t,Ro,Po|R') G,„(p)rfpdPo 
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It remains to uncouple the joint distribution of BP's path and momentum: 

l/o(i,Ro,Po|R') = Voit,AR;i^o) Git,Po\AR) , 
where last multiplier is conditional probability distribution of the momentum under 
given path value, and introduce conditional average correlation volume 

n„e<,(t,AR) = J y"r!„,g(t,AR,Po,p) G(t,Po|AR)G,„(p)dpdPo (33) 

After that inequality (|32p takes the form 

^0 + ^0 niaxF^'^^' AR) Vo{t, AR ; i^o) > (34) 

of declared restriction on rate of change of Vb(t, AR) . Before discussing it, consider 
another variant of the restriction [14] . 

3.2. Restriction on the BP's path distribution tails. Estimate 2 

Let n denotes at once a finite region in the r-space and volume of this region, 
presuming that it is centered near point Rq . Introduce n{S) = il{t, AR, Pq, p, S) as 
minimum of all those regions f2 what satisfy 



/ l/i(t,Ro,r,Po,p|R')dr - /l/i(i,Ro,r,Po,p|R')dr 
J yi(i,Ro,r,Po,p|R')dr 



< 



< 6 

with some fixed < (5 < 1 . That is flcori^) represents minimal region containing at 
least 100 (1 — S) percents of the total pair correlation. 

From the other hand, integrate identity (fT3|) over ri e il{S) . In view of the 
non-negativeness of Fi(f, Rq, ri, Pq, Pi |R') result must be non-negative: 

/ Ft'''\r\Ilo)dr K)(t, Ro, Po|R') G,„(p) + / Vi{t,Ro,v,Po,p\R.') dv > 

It is easy to verify that these two inequalities together imply a more interesting one: 

/ Fi(^^Hr|Ro)dr K)(i,Ro,Po|R')G™(p) + (35) 

J 0(5) 

+ {l-S) y"yi(t,Ro,r,Po,p|R')dr > 

Next, again let us replace here F^'^'^^ by max F^'^'''' (which even improves the 
inequality), then perform integration over momentums, apply relation (|25p and divide 
all by ( 1 - (5) . This yields 

1^0 max F^""^ n{t,AR, S) Vo{t,AR; i^o) + ?o ^^"ft'^^' > q (36) 

with characteristic average pair correlation volume defined by 

Tl{t,AR,S) ^ {1-S)-^ f f n{t,AR,Po,p,d) G{t,Po\AR)G,n{p)dpdPo (37) 



Unlike (I33p this characteristic volume has free parameter 0<(5<1. At (5^1 
almost all pair correlation is taken into account. But to make inequality (j36p most 
strong we have to minimize f2(t, AR, S) with respect to S . From this point if view, 
a choice when 1 — (S <C 1 can be preferred. At that, expectedly, ^{t, AR, S) is close 
to above defined Qneg{t, AR) if total pair correlation is negative. 
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3. 3. Finiteness of the correlation volume and failure of the Gaussian asymptotic 

Let us consider equilibrium and hence spherically symmetric Brownian motion what 
takes place in absence of the external force ( f = ). Assume that at sufficiently large 
scales, when t ^ t , with r being BP's mean free-ffight time, and 

(AR2(t)) = /AR2Vb(t,AR; i/o)dAR = 6Dt > SA^ , 

with A being BP's mean free path, the probability distribution of BP's path, 
Vo{t, AR ; uo) , tends to the Gaussian ([T]). Then it depends on gas density by way of 
the BP's diffusivity D only, and inequalities ([M]) or produce, together with (pS)) . 

where, respectively, 

ci(i, AR) = vq maxFi^'"'' f7„es(t, AR) , ci(t, AR) = vq max fI""^^ AR, 5) 

In case of gas, certainly, BP's diffusivity is a decreasing function of the density, 
c)lnI?/91nt/o < 0. Hence, the second addend in square bracket becomes negative 
at z = AR^/4£'i > 3/2 . As the consequence, inequality ([38|) can be satisfied if and 
only if at large values of z quantity ci (i, AR) grows at least proportionally to z . 

Notice that for not too dense gas (all the more, for dilute one) we can write 
A vtt and D k. vt^ = v'r''' 1 where vt \/T jM is characteristic thermal velocity 
of BP. Besides, for molecular-size BP (whose mass is comparable with atomic mass) 
Vt ^ Vs . This makes it obvious that limits of the quasi-uniform gas perturbation (see 
section 2.5) practically allows arbitrary large values of 2: , up to z ~ v^to/AD ^ to/r , 
where to is total duration of the "full-scale experiment". Therefore inequality ([55)1 
requires from ci{t, AR) and thus from i>qQ , with f2 being minimum of p3p and (j37p . 
ability to achieve as large values as to /t (where to , in its turn, is arbitrary large) . 

For dilute gas, more concretely, it is known that max F^'^''' = l ^ T Ovq/ dP — 1 , 
and A = (7rr^j/o)^^ oc \/vq, where rs is effective radius of BP-atom short-range 
repulsive interaction. Consequently, 91nZ3/91ni^o = —1 , and (f38|) satisfies only when 
1^0^^-1-3/2 > z = AB?/4:Dt at any z . 

In fact, we are faced with dilemma: either asymptotical statistics of 
Brownian path is Gaussian or volume of pair correlation is bounded above. 

What is better? In the first variant, "the law of large numbers" and conventional 
stochastic picture of Brownian motion hold true. But the strange enough requirement 
of unrestrictedly wide spatial statistical correlations indicates presence of self- 
contradictions in such theory. In the second variant, "the law of large numbers" 
does not work. But there are no nonphysical requirements and no contradictions. 

Undoubtedly, we have to prefer this second variant and claim that Brownian 
trajectory can not be imitated with the help of coin tossing or dice tossing or other 
"statistically independent" random events. 

3.4- Volume of pair correlation (volume of collisions) 

In order to make our consideration more pointed and a fortiori exclude from it any 
"collective" or "hydrodynamic" effects, a toy named "the Boltzmann-Grad limit" 
(BGL) is very useful. In this limit ^ 00 while rs ^ rA ^ {rs and rA 
are radii of short-range repulsive BP-atom and atom-atom interactions) in such way 
that gas non- ideality parameters 47rr^/3 and 47rr^/3 vanish but mean free paths of 
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of BP, A = (Trr^t'o)^^ , and atoms, A — (Trr^i^o)^^ , stay fixed. At that, BP collides 
with only infinitely small portion of atoms what surround it, therefore its previous 
collisions in no way can influence next ones. 

Folklore of kinetics includes opinion that at least under BGL the Boltzmann 
equation is true "zero-order approximation" for finite-density gas kinetics, and the 
corresponding Boltzmann-Lorentz equation (see e.g. [7]) for molecular BP (test atom 
or particle of rare impurity) is exact zero-order approximation for molecular Brownian 
motion. Since this equation inevitably yields [H [16] the Gaussian asymptotic, a best 
chance to resolve above formulated dilemma in detail is to consider pair correlations 
in the course of the BGL. For simplicity, again at f = . 

Recall a few things about pair CF, Vi(t, Rq, r, Pq, p|R') , already known from 
conventional theory [II O [8]. In respect to relative distance between two particles 
(in our case, BP and an atom), r — Rq , the pair correlation is accumulated inside 
a ''collision cylinder" which has radius « rs and is directed in parallel to relative 
velocity of the particles, v — Vq = p/m — Po/M . Importantly, characteristic value 
of the pair CF in this cylinder is comparable with product of one-particle DF. In 
our case this means that magnitude of Vi(t, Rq, r, Pq, p|R') (in particular, above 
defined quantity h{t, AR, Pq, p) ) is of order of first right-hand term in (fT3|) . As the 
consequence, magnitude of pair correlation, as well as Vo(<, Rq, Po|R') ; keeps safe 
under BGL, although inside more and more narrow collision cylinder only. 

But what is spread of pair correlation along the cylinder? Unfortunately, 
conventional theory never was interested in this issue, but in fact it reserves the spread 
to be infinite, at least at (v — Vq) • (r — Ro) > , i.e. for particles flying away after 
collision. As the consequence, integral / Vi{t, Ro, r, Pq, p|R') dr , integral in ([25l) and 
corresponding correlation volumes all turn to infinity. 

This became possible since the theory [SJ [8] neglected contributions from 
three-particle and other higher-order correlations, V2 , V3 , etc., although they involve 
collisions of the pair under attention with "third particles" (the rest of gas). That are 
not literally three-particle collisions but chains of (actual or virtual) pair collisions 
[TUl [TTl [TB] . It is not too hard to look after that magnitude of n-order CF on 
corresponding sets in n-particle phase spaces at n — 3,..., like at n = 2 , is 
of order of product of n one-particles DF irrespective to BGL. Therefore, due to 
collisions with "third particles" , the pair correlation of BP with atoms disappears when 
|r — RqI significantly exceeds A . Stronger separated particles hardly are participants 
of forthcoming or happening mutual collision, even if they aim precisely one to another. 

By neglecting all that, the theory unknowingly resolves the dilemma in favor of 
its first variant. That is why conventional kinetics, including the Boltzmann equation, 
is not true "zero-order approximation" and contradicts the virial relations. 

In fact, according to above remarks, a spread of the pair correlation along 
collision cylinder is finite and characterized by minimum of mean free paths A and A . 
Correspondingly, effective volume of the cylinder, or volume of pair correlation, O , is 
finite value of order of nr^A = v^^ (assuming, for simplicity, A ~ A), i.e. volume 
displayed per one atom. Since our reasonings are irrespective to where and when pair 
collisions take place, the estimate ~ Trr^A = Vq"^ is independent on t and AR 
and at once is estimate of upper boundary of the correlation volume. 

Let us consider this conclusion from the point of view of virial relations (fT5|) , 
with (221), (121, (ESI) and ^ with 
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3.5. God does not play dice with the Boltzmann-Grad gas 

It is easy to make sure that on the way to BGL, at any fixed function (/)(r) , expressions 
(HZl), (in]) and (USD simphfy to 

iy{r|0,Ro} v{vaA) , , , vo 
> 1 + 0(r) , > 1 + , > 1 , 

^{r,p|?/;,Ro} r-, , / ^ ^om 
> [1 + i/;(r,p)]G'„i(p) (39) 

Since change of vq during BGL is compensated by change of vb and va to 
keep constant A and A, left sides of (fT8|) . (f22|) . ([24]) and ([25]) become invariants of 
I'o although keeping their dependence on (j) . Consequently, all coeflicients in right- 
hand expansions into power series of (j) become independent on vq . On the left, 
the expansion means application of operators {d/d(j))"' which is equivalent to action 
of operators i'q [9/9j/o]" at fixed vb and va ■ Then under BGL any result of this 
operation also is independent on . 

What is for relations ^T^j and they turn to 

G^(p)Gm(Po) ^^^*^^'l^^'^o,Po} ^-f.[R'-Ro]/T ^ y,(t,R„,r,-Po,-p|R') (40) 
d^(r, p) 

with variational derivative taken at /i(r, p) — i'oGm(p). Obviously, this derivative 
represents reaction of probability distribution of BP's path, R' — Rq , to lodging, at 
t — , in point r one extra atom with momentum p . At that, importantly, it is 
presumed that this extra atom does not disturb initial equilibrium (of the statistical 
ensemble) , as if it was merely marked atom. 

Far enough under BGL, of course, this lodging can have an effect at such initial 
relative disposition of BP and the extra atom only which results in their direct collision. 
Hence, vector r — Rq should belong to those part of the collision cylinder which 
corresponds either to particles in m-state, i.e. flying one towards another, or to 
currently colliding (interacting) particles, i.e. separated by distance |r — Ro| < . 

Besides, |r — Ro| should not be much greater than min (A, A) . Otherwise 
collisions of either BP or extra atom with "third particles" (other atoms) will 
prevent the desired collision. This trivial notation once again, but more strikingly, 
highlights that volume occupied by the pair correlation on right-hand side of (|40)) 
can be estimated (at A ^ A) as fJ ~ Trr^A = Vq"^ irrespective to the BP's path. 
Indeed, on the left in ()40|) . in view of the causality principle, a whole future path 
of BP can not influence possibility of its collision with concrete marked atom at 
very beginning of the path. Moreover, the mentioned path and collision can not 
be mutually statistically correlated since the marked atom and BP deal with different 
non-intercrossing collections of "third particles". Therefore, any dependence of the 
right-hand pair CP on t and AR characterizes its magnitude but not its spread. 

Now we can say almost with confidence, that correlation volumes 0„eg(i, AR) , 
ri(t, AR, Pq, p, (5) = ri((5) and f7(t, AR, (5) , introduced in sections 3.1 and 3.2, have 
upper boundary ^ v^^ independent on t and AR or z , and should be written 
simply as l^neg , ^^(Po,P, (5) = ^{5) and ^{5). Correspondingly, ci(t, AR) in ((38)) 
is not a function but a constant of order of unit. 

Thus, the second resolution of the dilemma from section 3.3 is acceptable only. We 
can say that the god of mechanics does not play dice. He disposes of stronger means: 
Hamiltonian dynamics of many particles {N > to/r , according to [T31 HH [H] ) is able 
to create much more rich randomness than dice tossing can do. 
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3.6. Probability distribution of BP's path possesses power-law long tails 

Now we are ready to discuss what kind of asymptotic of BP's path distribution is 
really allowed instead of the Gaussian one. 

According to previous section, let us rewrite inequahties l|34|) or ((36|) in the form 



c,Vo{t,AR;.o)+T^o ^^'^\^^''''''^ >0, (41) 
Here ci is a constant whose upper estimates obtained in two different ways look as 
ci = vq ma^F^'^'^^ ^neg , c-i = max f['^'^'^ Q,{5) 

with ^Ineg and r2((5) being volumes of pair correlation (volumes of collision) defined 
by (|33p and (I37p and bounded above by a value ~ i/q^ (volume per one atom), so 
that ci ~ 1 . Probably, there are methods to estimate ci differently from sections 3.1 
and 3.2. Anyway, to make inequality (|4ip stronger, we should take minimum of all 
available estimates. 

For simplicity, consider equilibrium Brownian motion at f = . In case of 
dilute enough gas, or even liquid in three or more dimensions, when hydrodynamical 
contributions to BP's velocity are far from domination, and the "diffusion law" 
AR^ (X t holds, it seems natural if asymptotically, at i ^ r , the BP's path 
distribution Vb(t, AR;i'o) is characterized, similarly to VG'(i, AR), by a single 
parameter, that is diffusivity. For 3-D space, 

Vo{t, AR ; vq) (4L>t)-3/2 q>{AR^/4Dt) (42) 
Here it is presumed that 

J *(a2)da = 1 , J a2*(a2)da =3/2 (43) 

The first of these requirements is the normalization condition, while the second means 
that (AR2(t)) = 6Dt and in fact serves as quantitative definition of the diffusivity. 
Then, inequality (|4T|) yields 

- 1 



vo dD 

if diffusivity is a decreasing function of density, dD / dv^ < , and 



d*(z) 3 
a^(z) + z > , a=-+ci 



(44) 



a^>(z) + z — A-L < , a = a 

dz 2 



vq dD 
15 



■ 1 



in opposite case dD/du^ < . 

Confining ourselves by the first case, we see that function ^'(z) must have power- 
law long tail: '^{z) (x z~^ at z — )■ cx3 , where (3 < a . At that, in order to satisfy 
second of conditions ((43l) . exponent a should be greater than 5/2 . 

Example of such behavior, with maximally possible (3 — a , is presented by 

= , y„(t,AR;.o)- -^^fl + ^V", (45) 

(1 + ^)" ' ' ' (4D'i)3/2 1^ AD'tJ ' ^ ' 

where 7(a) = 7r-3/2 r(c,)/r(a - 3/2) and D' ={a-5/2)D. 

Unpleasant aspect of such behavior is unboundedness of high enough statistical 
moments of AR . However, recall that real Brownian motion in addition to diffusivity 
has at least one more important parameter, namely, the BP's thermal velocity vt ■ 
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Then Vo{t, AR) should be considered as a function of two dimensionless parameters, 
z — AR^/iDt and y = "R? /v'^i^ , and asymptotically instead of (M)) we can write 

Here ^'(z) as before satisfies while function Q{y) w 1 at y ^ 1 and fast enough 
tends to zero at ?/ > 1 thus sharply cutting Vb(i, AR) when |R| > vrt . This variant 
also is allowed by inequality (|4ip . because vt is independent on density. And now all 
statistical moments of BP's are finite. 

Nevertheless, from the point of view of higher-order moments, such asymptotic 
still is not exhaustive, excepting case of dilute gas. For sufficiently dense gas (all the 
more, liquid) we should take into account, at the minimum, such third (density- 
dependent) parameter as isothermal speed of sound, = yTvo/ravQ , with vq 
presented by (p6|) . At the same time, such complications practically do not change 
main probabilistic characteristics of BP's path, as well as its mean square behavior, 
and in no case cancel the long power-law tails. 



3. 7. Comparison with solutions of the BBGKY equations 



Previous results, obtained from the first-order virial relations (one- and two-particle 
CF) only, are in remarkable qualitative (and even semi-quantitative) agreement with 
results obtained in [TB] from infinite chain of roughened BBGKY equations describing 
a test (marked) gas atom in the role of BP in the Boltzmann-Grad gas. But, of 
course, our approach here could not give numeric value of exponent a in (|44p and 
(|45p . Approximate analysis in [T^] gave distribution ((21), that is a = 7/2, which 
corresponds to ci — = 2 . 

More precisely, solution obtained in [16] indeed has the form (|46| . At that, the 
cut-off function Q works already at fourth-order statistical moment, yielding 

( AR4(i) ) = j AR4 Vo{t, AR) dAR ^ 3 ( AB?{t) ? (47) 



with (AR2(t)) = &Dt. 
quantity, D , with mean 



This result looks as if BP's diffusivity was fluctuating 
{b) = D and variance {&) - (D)'^ w \n{t/T) which 
depends on total duration of BP observation, t . 

From the point of view of ^ variance of the fluctuating diffusivity, D , is infinite. 
But probability distribution of D is quite certain. One uncovers it representing ^ 
as superposition of Gaussians with various values of diffusivity: 

- 7/2 



Vo(i,AR) 



r(7/2) 

{4T:Dt)y^ 
wiC) 



1 



AR^ 



4Dt 



1 



{AT:CDt)y^ 



exp 



AR 



2 1 



ACDt 



exp 



(48) 



where ^ represents D/D while w(C) is probability density of C- 

At one and the same D distribution ^ or (j48|) seems rather dissimilar to the 
Gaussian (jT]) . But in fact the only essential difference between them is long tail of ^ . 
To make this visible, we may rescale ([2]) by suitable increase ( "renormalization" ) of 
its diffusivity and compare one-dimensional projections of ([T]) and ^ onto some axis 
X . Left plot on Figure 1 illustrates such comparison, with rescaling factor Sy/ir/i 
which equalizes heights of both distributions at X = . 
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2 4 6 2 4 

X = X/(2Dt)'^2 x=XT/Dft 

Figure 1. (a) Projections onto X-axis of reduced Gaussian distribution JTJ, 
(27r)^^/^ exp (— (thin black curve), and reduced rescaled non-Gaussian 
distribution l[2]l, {2tt)~^^^ {l + 8x^ /9tt)~^^^ (thick gray curve); the inset magnifies 
tails of the distributions, (b) Reduced distribution (|51^ . x~^exp(— 1/x) (thick 
gray curve), of drift displacement in comparison with reduced rescaled result, 
5{x — 16/97r) (black vertical line), of usual theory. 



Although in [16j only thermodynamically equilibrium Brownian motion was 
considered, the integral representation as combined with FDR (PO)) . can be used 
to predict statistics of non-equilibrium Brownian motion at f 7^ . Reasonableness 
of such operation, for a weak non-equilibrium, was confirmed e.g. in |2H I23j. The 
word "weak" means that BP's drift velocity is much smaller than vt and and 
therefore is linear function of the external force: 

(AR(t)) = / AR Vb(t, AR)dAR = Vdt , = Dt/T 

Here D/T is mobility of BP, in agreement with the Einstein relation [HIS which, by 
the way, directly follows from ^ or ^ (see also [TTlIinilll])- 

To underline dependence of Vb(t, AR) on the external force, let us denote it 
as V/(t, AR) . Then, if we do not wish to look on so far distances as vxt (all the 
more, |f|t^/2M), we can find asymptotic of V/(t, AR) at t ^ r in the form (|48|) 
but replacing equilibrium Gaussian components by non-equilibrium ones: 



WAR)-/^^^^exp 



[AR-CDtt/T] 



w{0 dC (49) 



Recollecting that d\nD/d\ni^Q = —1 under BGL, it is easy to make sure that this 
expression satisfies inequality (|4T|) , that is 



c,V,it,AR)-D^M:^>0 , = 2 , 

as well as equality (|20|) . In short, (|49|) agrees with both virial relations and FDR. 

Consider (|49ll at a time when drift of BP exceeds its diffusion, |vd|t > V 6Dt . In 
other words, when f -v^ t ^ T , i.e. work of the force much exceeds thermal energy per 
degree of freedom. At that, distribution Vf{t,A'R) becomes highly asymmetric and 
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anisotropic. Therefore we have to distinguish BP's displacement X afong direction 
of the force f and two orthogonal displacements Y and Z . In such notations 
AIL = {X, Y, Z} , and ^ yields 



Vf{t,{X,Y,Z}) ^^^exp 



|f|(r2 + z2 



'd 



r { |Vd|t 



X3 "^Kx 



(50) 



ATX 

At any fixed Y and Z this expression has long tail in positive X-direction, so that 

Vf{t,X) ^ j j Vj{t,{X,Y,Z})dYdZ = ^ exp (-^) , (51) 

but in strictly this direction only. Right-hand part of Figure 1 shows this asymptotical 
distribution of BP's drift, in comparison with drift distribution in conventional 
kinetics, which is nothing but delta-function S{X — \'Vd\t) ■ At the same time, 

JVf{t,{X,Y,Z})dX = JVQ{t,{X,Y,Z})dX = {2TTDt)-^ [1 + {Y^ + Z^)/'iDt]-^ 
has long tail in any direction in YZ-plane, and anybody who keeps under observation 
Y and Z only sees no signs of BP's drift along X-axis. 

Basic idea of [TQ, earlier suggested in [l0l[TT]), was to consider DF averaged over 
"collision boxes". The latter are just sets in n-particle phase spaces already mentioned 
in section 3.4. For n = 2 , "collision box" is nothing but above exploited "collision 
cylinder" . In general, "collision box" is a "skeleton" set of n-particle configurations 
snapped up from chains of n — 1 connected (actual or virtual) collisions, under given 
input momenta. Thorough description of all these sets hardly is possible, but hardly 
it is obligatory. Much more important thing is that even under rough description one 
immediately discovers that the Boltzmann-Grad gas is true terra incognita . 

3. 8. Origin of the historical correlations and long tails and 1 // noise 



Let us return to identity (|40|l . One and the same vector r — Rq on its opposite sides 
belongs to opposite half of the collision cylinder. Hence, pair correlations on right- 
hand side are concentrated at out-st&tes (i.e. concern particles flying away one from 
another) and at central region of the cylinder (i.e. concern also currently interacting 
particles, with |r — Ro| < ), while m-states are uncorrelated. 

This consequence of (|40|) is in agreement with usual theory [5j [6l [71 [8] . The latter 
considers it as a good reason to believe in "statistical independency" of colliding 
particles (Boltzmann's "molecular chaos"). In present theory, however, actually 
colliding particles occupy central part of collision cylinder where, according to ([^0]) . 
Vi{t, Ro, r, Pq, p|R') is non-zero, and thus these particles are mutually correlated and 
statistically dependent. Hence, statistical dependency is not "exported from outside" 
but occurs at the time of collision. 

The point is that Vi(i, Rq, r, Pq, p|R') (as well as higher-order CF) represents 
"historical correlations" (see section 2.3): it treats current collision as not a separate 
event but last term of long random sequence of BP's collisions with atoms on its way 
from R' to Rq . Therefore statistical dependencies involved by Vi(t, Rq, r, Pq, p|R') 
should be addressed not to currently colliding particles (BP and atom) themselves but 
to long-range fluctuations in relative frequency of BP's collisions. 

With the purpose to make certain that this is only reasonable understanding of 
the correlations, in particular spatial ones, Vn{t, AR \ vq) ( n > 1 ), let us introduce, 
in terms of section 3.2, function 



W^i(t,AR) = 



/ i^i(t,Ro,r,Po,p|R')rfr 

JO(Po,p,5) 



dpdPo 
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which is direct analogue of function W2 from [TB] and satisfies 

J Wi{t, AR) dAR = ci lyoUiS) 

Evidently, AR) is probability density of finding BP at point AR (after start 

from coordinate origin) and simultaneously some atom in its close vicinity in some 
m-state or owi-state of their mutual collision. Consequently, ratio 

Wi(i,AR)/Fo(t,AR) = Vpostit^AK) 

presents a measure of conditional probability of BP's collision under condition that 
path AR is known. The subscript "post" underlines that in essence this is a posteriori 
probability. In absence of historical correlations between BP and atoms we would have 
Wi(t, AR) = ciVo(t, AR) and above ratio would reduce to unconditional a priori 
probability measure: Vprior = Vt{5) . 

In conventional kinetics, based on concepts like "probability of collision", any 
difference between Vpost and Vprior is unthinkable since it means dependence of 
"probability of collision" on wherefrom BP had started in former times. Therefore 
factual difference of Vpost from Vprior , indicated by the virial relations, says that a 
random point process formed by BP' collisions in not a usual Poisson process. Then, 
what is it? The answer is rather obvious: it is Poisson process but supplied with 
such (relatively slow) scaleless random variations of "probability of collision" which 
exclude possibility to get it by time averaging. Indeed, on second thought it is clear 
that Poissonian statistics introduces "too strict disorder" in collisions to be likely. 

Such interpretation of inter-particle correlations which keep safe even under BGL 
was expoimded in [lOl [11] (besides, it was in part developed in [131 [H [HI [16] . 
and its most principal aspects were anticipated already in [2Ql [211 [221 [23]). It was 
demonstrated that fluctuations in "probability of collision" (or, more precisely, relative 
frequency of collisions) possess scaleless 1/f-type spectrum (i.e. represent "1/f-noise") 
and may be described also as 1/f fluctuations of BP's diffusivity (and mobility). Now 
we arrived at principally same results after start from virial relations. 

Scaleless character of the fluctuations is due to fact that the system under 
consideration constantly forgets history of BP's collisions and therefore has neither 
stimulus nor means to force relative frequency of collisions to be certain. To keep it 
certain at arbitrary large time scales, the system, in opposite, should keep in mind 
arbitrary old history. In this sense, a priori equalization of Vpost and Vprior acts 
as infinitely long memory. Notice that in comparison with Poissonian statistics of 
collisions a real one to some extent resembles Bose statistics. 

Now let us compare Vpost and Vprior ■ From section 3.5 it follows that sign of 
Vi(t, Ro, r, Po, p|R') is determined by t and AR only. Taking this into account, in 
the spirit of section 3.2 we can derive one more inequality. 



W^i(i,AR)-ciFo(i,AR)- ^i(*'^^) 



< 5 



Vi{t,AR) 



1 



(52) 



l-S 

Being combined with ([25| it implies that 

Vpostjt^AR) < ^ ^ 1 d\nV^{t,AR) ^ ^ 

1^ prior 

ci d In vq 

when aVo(i, AR)/at'o < . When, in opposite, SVb (t, AR)/at/o > , then both 
inequality signs in (I53p should be inverted. 

Inequality ([55)1 represents strengthened form of (|25p . At tails of Vo(t, AR) always 
9Vo(i, AR)/9fo < 0, and it shows that, naturally, a posteriori probability of BP's 
collisions far enough at tails is smaller than a priori one. 
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Notice that if ci was an exact value, determined by minimum of all possible 
estimates of pair correlation volume (see section 3.6 and above), then middle expression 
in (j53p would achieve its theoretical minimum too. Hence, in fact left inequality should 
be replaced by equality. Similar conclusion is valid also in respect to its alternative at 
dVo{t,AR)/di^o > 0. Thus we can write 

Vpostjt, AR) ^ ^ ^ 1 d\nVo{t,AR) ^ ^ _ 1_ ginFoft,AR) 

1^ prior 

The last equality here relates to the Boltzmann-Grad gas. For distribution ([2]) or ([48]), 



Vpo,t{t, AR) ^ 7 / ^ AR^Y 

1-^ prior 4 y ADt J 



(55) 

prior 

So substantial changeability of the a posteriori "probability of collision" shows that 
in fact this quantity has no certain value obtainable by time averaging. 



4. Conclusion 



Being guided by wrong idea of "statistical independency" of colliding particles 
(Boltzmann's "molecular chaos"), classical gas kinetics neglected statistical inter- 
particle correlations which inevitably arise in spatially non-uniform Gibbsian 
statistical ensembles. Thus the concept of a priori definable "probability of collision" 
was unknowingly imposed on the theory as characteristics of a concrete particle 
trajectory. As for Brownian motion (self-diffusion) of gas particles, the result was 
"the law of large numbers" stating that probability distribution Vb(i, AR) of path of 
molecular Brownian particle (BP) is drawn towards the Gaussian distribution. 

We derived exact virial expansion connecting response of Vb(<, AR) to gas 
perturbations, from one hand, and pair and many-particle "historical" statistical 
correlations between the BP's path and gas, from the other hand. Specificity of 
"historical" correlations is that they are just products of initial spatial non-uniformity 
of an ensemble. Thus existence of such correlations is proved. At the same time, 
we showed that finiteness of spatial spread of these correlations is incompatible with 
Gaussian asymptotic of Vb(i, AR) even in dilute gas (under the Boltzmann-Grad 
limit). Then we demonstrated that the spread is really finite and, consequently, 
Vb(iiAR) has essentially non-Gaussian form, with power-law long tails (lasting up 
|AR| ~ vxt). These results mean (i) that BP's path can not be divided into 
"statistically independent" events or pieces, and, as combined with our previously 
published [M UHl HI 121] and unpublished [TTl [H [H [TH [H] results, (ii) that 
"probability of collision" or, equivalently, diffusivity (and mobility) of BP in fact 
undergo scaleless fluctuations like the so-called 1/f-noise. 

In principle, these conclusions extend to the kinetics as the whole. If applied 
to charge carriers in semiconductors, they naturally explain the "inherent 1/f-noise" 
under experimental investigation during many years [251 126] . However, one has to 
revise kinetics of electron-phonon systems (as well as phonon systems themselves 
[27]), again removing from it a priori "probability of collision" and "statistical 
independencies" , and returning to honest analysis of infinite chains of many-particle 
statistical correlations. Perhaps, a time of such revision is not far off. 

I am grateful to Dr. Yu. Medvedev and Dr. I. Krasnyuk for useful discussions. 
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